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Abstract 



Within the effective Lagrangian framework, we explicitly evaluate the partition 
function of two-dimensional ideal ferromagnets up to three loops at low temperatures 
and in the presence of a weak external magnetic field. The low-temperature series 
for the free energy density, energy density, heat capacity, entropy density, as well 
as the magnetization are given and their range of validity is critically examined in 
view of the Mermin- Wagner theorem. The calculation involves the renormalization 
and numerical evaluation of a particular three-loop graph which is discussed in detail. 
Interestingly, in the low-temperature series for the two-dimensional ideal ferromagnet, 
the spin-wave interaction manifests itself in the form of logarithmic terms. In the free 
energy density the leading such term is of order T 4 In T - remarkably, in the case of the 
three-dimensional ideal ferromagnet no logarithmic terms arise in the low-temperature 
series. While the present study demonstrates that it is straightforward to consider 
effects up to three-loop order in the effective field theory framework, this precision 
seems to be far beyond the reach of microscopic methods such as modified spin-wave 
theory. 



1 Introduction 



In a very recent article on the thermodynamic properties of two-dimensional ideal 
ferromagnets [lj], the general structure of the low-temperature series for the free energy 
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density has been discussed using the method of effective Lagrangians. In particular, 
it has been argued that the spin-wave interaction does not yet manifest itself at 
order T 3 in the free energy density - it only shows up at order T 4 lnT. The explicit 
evaluation of the various Feynman graphs contributing at this order, however, has 
not been performed in that reference - this will be the subject of the present article. 
The calculation involves a particular three-loop graph whose renormalization and 
subsequent numerical evaluation - although quite elaborate - nevertheless is rather 
straightforward within the effective field theory framework. While Ref. [lj focused 
on the low-temperature series for the free energy density, here we also consider the 
magnetization and critically examine the range of validity of the corresponding low- 
temperature series - indeed, in view of the Mermin- Wagner theorem, one has to be 
very careful by taking the limit of a zero external magnetic field H. 

Interestingly, these logarithmic contributions T n In T in the partition function are 
restricted to two-dimensional ferromagnets and do not show up in the case of three- 
dimensional ideal ferromagnets. There the low-temperature series of the various ther- 
modynamic quantities consist of integer and half-integer powers of the temperature 
Actually, the occurrence of such logarithmic contributions is well-known in the 
domain of particle physics - the low-temperature expansion of e.g. the free energy 
density in Quantum Chromo dynamics also involves a logarithmic term T 8 lnT j^j. A 
term T 8 In T also occurs in the low-temperature expansion of the free energy density 
of the three-dimensional antiferromagnet Q, which obeys a linear, i.e. relativistic, 
dispersion law. 

We also confront our series obtained within the effective Lagrangian framework 
with the relevant literature and point out that our approach is much more efficient 
than conventional condensed matter methods such as spin-wave theory or Schwinger- 
Boson mean field theory. In particular, while all these different studies were restricted 
to the idealized picture of the free magnon gas, in the present article we discuss in 
detail how the spin-wave interaction manifests itself in the low-temperature behavior 
of the two-dimensional ideal ferromagnet. 

The rest of the paper is organized as follows. In Sec. [2]we briefly discuss some essen- 
tial aspects of the effective Lagrangian method - detailed accounts can be found in the 
references provided below. The evaluation of the partition function up to three-loop 
order in the low-temperature expansion is presented in Sec. [3j While the renormaliza- 
tion up to order p 6 oc T 3 turns out to be straightforward, the handling of ultraviolet 
divergences at order p 8 is more involved and is considered in detail in Sec. HI The 
low-temperature expansions for the free energy density, energy density, heat capacity, 
and entropy density of the two-dimensional ideal ferromagnet are given in Sec. El The 
low-temperature series for the magnetization is discussed in Sec. [6] and the range of 
validity of this series is critically examined in view of the Mermin- Wagner theorem. 
Here we also compare our results with the condensed matter literature. While our 
conclusions are presented in Sec. the numerical evaluation of a specific three-loop 
graph, the logarithmic renormalization of effective constants and the representation 
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for the pressure at zero magnetic field are relegated to three separate appendices. 



Unfortunately, the systematic and model-independent effective Lagrangian method 
is still not very well known within the condensed matter community. Therefore we 
would like to provide the reader with a brief list of references where condensed matter 
problems have been successfully solved within the effective field theory framework. 
These include antiferromagnets and ferromagnets in two and three spatial dimensions 
, as well as two-dimensional antiferromagnets which are the precursors of 



high-temperature superconductors [20-29 



Many different tests - both analytical and numerical studies - unambiguously 
confirm that the effective Lagrangian technique represents a rigorous and systematic 
framework to address condensed matter systems which exhibit a spontaneously broken 
symmetry. On the one hand, the correctness of the effective Lagrangian approach 
was demonstrated explicitly in Ref. [3o| by comparing the microscopic results of an 
analytically solvable model for a hole-doped ferromagnet in 1+1 dimensions with 
the effective theory predictions. On the other hand, in a series of high-accuracy 
investigations of the antiferromagnetic spin- j; quantum Heisenberg model on a square 



lattice using the loop-cluster algorithm [3 lM34| , the Monte Carlo data were confronted 



with the analytic predictions of the effective Lagrangian theory and the low-energy 
constants were extracted with permille accuracy. 



2 Essential Aspects of Effective Lagrangians 

The present article deals with the explicit evaluation of the partition function of the 
two-dimensional ideal ferromagnet at three-loop order and the subsequent discussion 
of various thermodynamic quantities - including the magnetization - at low tem- 
peratures. It is an extension of the work presented in the recent article [H, where 
the evaluation was restricted to two-loop order. The basic principles of the effective 
Lagrangian method and the perturbative evaluation of the partition function have 
been outlined in section 2 of that reference and will not be repeated here in detail. 
The interested reader may also find a detailed account on finite-temperature effective 
Lagrangians in appendix A of Ref. j^] and in the various references given therein. In 
addition, for pedagogic introductions to the effective Lagrangian technique, we refer 



to Refs. [35H39|. In this section, we will just focus on some essential aspects of the 



effective Lagrangian method at finite temperature. 

The effective Lagrangian, or more precisely, the effective action 

S eff = J d 3 xC eff (2.1) 

of the two-dimensional ideal ferromagnet must share all the symmetries of the under- 
lying Heisenberg model, i.e. the spontaneously broken spin rotation symmetry 0(3), 
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parity and time reversal. The basic degrees of freedom of the effective Lagrangian are 
the two real magnon fields - or the physical magnon particle - which emerge due to 
the spontaneously broken spin symmetry 0(3) — > 0(2). 

The various pieces in the effective Lagrangian are organized according to the num- 
ber of space and time derivatives which act on the magnon fields. This derivative 
expansion - or expansion in powers of momentum - is completely systematic: At low 
energies, terms in the effective Lagrangian which contain only a few derivatives are 



the dominant ones, while terms with more derivatives are suppressed |40H42| . The 
leading-order effective Lagrangian for the ideal ferromagnet is of momentum order p 2 



and takes the form 11 



d ff = tf'f ™ + ^HU 3 - l -F 2 d r U*d r U l . (2.2) 

The two real components of the magnon field, U a (a = 1, 2) are the first two compo- 
nents of the three-dimensional unit vector U l = (U a , U 3 ). The quantity H is the third 
component of the magnetic field H = (0, 0, H). While the derivative structure of the 
above terms is unambiguously determined by the symmetries of the underlying the- 
ory, the two a priori unknown low-energy constants - the spontaneous magnetization 
at zero temperature S, and the constant F - have to be determined by experiment, 
numerical simulation or comparison with the microscopic theory. Note that one time 
derivative (d ) is on the same footing as two space derivatives (d r d r ), i.e., two powers 
of momentum count as only one power of energy or temperature: k 2 oc u, T. 

The next-to-leading order Lagrangian for the ideal ferromagnet is of order p 4 and 
amounts to 0] 

C 4 eff = hidrWdrU'f + hidrWdsU*) 2 + l 3 AW AW , (2.3) 

where A denotes the Laplace operator in two spatial dimensions, d s — 2, and the 
quantities lx,l% and 1% are effective coupling constants. 

A crucial point underlying the perturbative evaluation of the partition function 
concerns the suppression of loop diagrams in the effective field theory framework. 
In two spatial dimensions, ferromagnetic loops are suppressed by two powers of mo- 
mentum. This suppression rule lies at the heart of the organization of the Feyn- 
man graphs of the partition function for the two-dimensional ferromagnet depicted in 
Fig. [TJ For example, the two-loop diagram 8d with an insertion from £. 4 ff, containing 
four magnon fields, is of order p 8 , as it involves (p 4 ) and two loops (p 4 ). 

An inspection of the diagrams in Fig. [T] reveals that insertions from only 
appear in the one-loop graph 8g. The only term we need for the present evaluation is 
thus quadratic in the magnon field, 

C 6 eff = alPA 3 ^ , (2.4) 



4 




8/ 8g 



Figure 1: Feynman graphs related to the low-temperature expansion of the partition 
function for a two-dimensional ferro magnet up to order p 8 . The numbers attached to 
the vertices refer to the piece of the effective Lagrangian they come from. Vertices 
associated with the leading term C? e ^ are denoted by a dot. Note that ferromagnetic 
loops are suppressed by two powers of momentum in two spatial dimensions, d s =2. 
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where C\ is an additional effective low-energy constant. 



Lattice anisotropics do not yet manifest themselves in the leading-order effective 
Lagrangian C? e ^ - the expression (12.21) is invariant under continuous space rotations, 
even though the underlying lattices (square, honeycomb, triangular and Kagome, to 
name the most prominent ones) are only invariant under discrete space rotations. 
Although this accidental symmetry is restricted to £e//> m the present analysis we 
assume that the higher-order pieces and are 0(3)-space-rotation invariant 
as well. The conclusions of the present paper will not be affected by this idealization. 



3 Partition Function up to Three Loops 

While the evaluation of the partition function of the two-dimensional ideal ferromag- 
net up to order p & was presented in Ref . [l| , here we go one step beyond and consider 
the evaluation up to order p 8 , where three-loop graphs start to show up. This calcu- 
lation is rather involved due to the renormalization and numerical evaluation of the 
three-loop graph 8c. 

As depicted in Fig. [H there are only four Feynman diagrams that contribute to 
the partition function up to order p 6 or, equivalently, up to order T 3 . The evaluation 
of these diagrams can be found in Ref. [l| - the final expression for the free energy 
density amounts to 

^-^-^r^—-^ T ^— + o(v&) < (3 - i) 

71=1 71=1 

where (3 = 1/T and H is the external magnetic field. The terms of order T 2 and T 3 
arise from one-loop graphs and thus represent contributions to the free energy density 
originating from noninteracting spin waves. While the former term is exclusively 
determined by the leading-order effective constants £ and F (7 = F 2 /E), the term 
of order T 3 , on the other hand, involves the next-to- leading-order effective constant 
Z 3 from jC-eff- What is quite remarkable is that the spin- wave interaction does not 
yet manifest itself at this order in the low-temperature expansion of the free energy 
density. The only potential candidate, the two-loop diagram 6a of order T 3 , turns out 
to be zero due to parity [l|. 

Before we evaluate the partition function at order p 8 , we have to recall that, at 
finite temperature, the propagator is given by 

00 

G{x) = A(f , x 4 + , (3.2) 

n = —00 

where A(x) is the Euclidean propagator of ferromagnetic magnons at zero tempera- 
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ture, 



A(x) 



dk A d 2 k e ikx-ik iXi 



(3.3) 



(2tt) 3 7 ^2 _ iki + ■ 

An explicit representation for the thermal propagator, dimensionally regularized in 
the spatial dimension d s , is 



1 d s OO 1 



(2tt)^ V 7 , — ^ 

v ' n = — oo In 



cxp 



ar 



fxHx n 



with 



x n = x 4 + n/3 . 



e(x n ) , (3.4) 



(3.5) 



Also, it is convenient to introduce the following notation, 



G(x) 



x=0 



G A 



AG(x) 



G An 



X=0 



A n G(x) 



x=0 



(3.6) 



where A represents the Laplace operator in the spatial dimensions - note that the 
quantity A(x), on the other hand, stands for the zero-temperature propagator. 

The quantities G±, Ga, as well as thermal propagators involving even higher-order 
space derivatives, can be split into a finite piece, which is temperature dependent, 
and a divergent piece, which is temperature independent, 



Lri — Lr 1 + Lr 1 , Lr A — l^A + ■ 

The explicit dimensionally regularized expressions are 

1 



(3.7) 



and 



Gl = 



G° A 



Gl = 
G\ = 

1 

(2tt)^ 
1 

(2tt)^ 



(27r) d « 
1 

(2vr)^ 



7 



E 



7T\ 2 



7 



n=l M) 
1 

^7 ex P 

X 4 2 



2 



->2 

ar 



47X 4 



8(2:4) 



(3.8) 



x=0 



7 



7T\ T" 
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^ ^ (n/3)f + 1 

I I — x 2 

+ 



x fi +1 1 2 7 ' 4 7 % 4 



cxp 



47x4 



e(x 4 ) 



(3.9) 



x=0 



The temperature-independent pieces G\, G A , . . . are all related to momentum inte- 
grals of the form 



J d d °k(k 2 y 



cxp 



— 7X4A; — x^jiH , m = 0,1,2,..., 



(3.10) 



which are proportional to 



exp[-x A fiH] ( d s 
— 1 [m H 

( 7 x 4 ) m+ ^ V 2 



(3.11) 



In dimensional regularization these expressions vanish altogether: Gi,G° A , and zero- 
temperature propagators involving higher-order space derivatives do not contribute 
in the limit d s — > 2. The only contributions we are thus left with here, are those which 
involve the temperature-dependent pieces. 

According to Fig. [TJ there are seven diagrams at order p 8 we need to evaluate. We 
first consider the two one-loop graphs which contain vertices from £<!// an d ^e//- For 
graph 8g, which only involves an insertion from £e//' we obtain 

2d 



(3.12) 



yielding the temperature-dependent contribution 

rp 3 Cl 



r 4 E 



^ e -nHn/3 



n=l 



n 1 



(3.13) 



Graph 8f, which contains two insertions from £e/p * s proportional to an integral over 
the torus T = 1Z da x S 1 , with circle S 1 defined by —f3/2 < x 4 < (3/2, and involves a 
product of two thermal propagators, 



z sf 



2l\ 



Jd ds+1 xA 2 G(x) A 2 G(-x). 



(3.14) 



Integrals of this type, as shown in Ref. [4J, can be converted into an expression dis- 
playing one propagator only, 



- I d ds+1 y A rn G(-y) A n G(y) 



A(m+n) dG ( X ) 



(3.15) 



x=0 



such that one ends up with 



Z 8f 



211 



A 



4 dG(x) 
d(fiH) 



(3.16) 



x=0 



Accordingly, the temperature-dependent part of graph 8f amounts to 



12/f 



4 8/ 



7rS 2 7 5 



r 4 E 



x> e ~/iHn/3 



n=l 



(3.17) 



We now turn to the two- loop graphs which involve insertions from anc ^ ^"eff- 
Graph 8d contributes with 



2 21 
z sd = ~^( 8/ i + 6/ 2 + 5/ 3 ) G A G A - G 1 G A 2 . 



(3.18) 
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The evaluation of graph 8e yields 



Z8e = -^f GiGa2 ' (3.19) 

and cancels the second term in zgd- For the temperature-dependent part originating 
from the two-loop graphs of order p s we thus end up with 



z. 



8h + 6l 2 + 5l 3 m4 Kp, e 



T 

8[de] — 0/1„2V2o,4 



N -nHn/3 " 2 



24tt 2 S 2 7 4 ^ n 

1 n=l 



Finally, we consider the three-loop graphs. Note that they exclusively contain 
vertices from the leading-order Lagrangian C? e ^- Graph 8a factorizes into a product 
of three thermal propagators (and space derivatives thereof), to be evaluated at the 
origin, 

z 8a = -|!g a (G 1 ) 2 . (3.21) 
Remarkably, the three-loop graph 8b turns out to be zero, 

z sb = . (3.22) 

Much like the two-loop graph 6a, the three-loop graph 8b does not contribute to the 
partition function. 

We are left with the cateye graph 8c which leads to 

z 8c = ~K + ^GaCGO 2 . (3.23) 

The quantity K denotes the following integral over the torus involving a product of 
four thermal propagators 

K = Jd ds+1 xd r Gd r Gd s Gd s G. (3.24) 

Here we have used the notation 

G = G{x) , G = G(-x) . (3.25) 

The second term in (I3.23P cancels the contribution from graph 8a - hence the overall 
contribution from the three-loop graphs is just the one proportional to the integral K. 
Note that, unlike all other pieces in the free energy density up to order p 8 , this quantity 
is not a product of thermal propagators (or derivatives thereof) to be evaluated at 
the origin - its structure is much more complicated. In the next section, as well as in 
the appendices \M and [BJ we will address in detail the renormalization and numerical 
evaluation of this integral which contains a total of four infinite sums. 
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4 Renormalization of the Cateye Graph 



In order to renormalize the three-loop graph 8c, we will follow the method outlined 
in Ref. [f|, where the same graph was considered in a Lorentz- invariant framework. 

The relevant integral from the three-loop contribution is given by 

K = Jd ds+1 xd r Gd r Gd s Gd s G, 

and involves a product of four thermal propagators. In order to analyze the limit 
d s — >2, the thermal propagators are split into two pieces, 



G(x) = G T {x) + A(x) . 



(4.1) 



While the zero-temperature propagator A(x) contains the ultraviolet singularities, 
the temperature-dependent part G T (x) is finite as d s — >2. Note that, if we restrict 
ourselves to the origin, we reproduce the first relation of Eq. fl3.7p . 

With the above decomposition, the integral K yields nine terms that can be 
grouped into the following six classes - for simplicity we do not display the derivatives: 



A : 


G T (x) G T {x) G T {-x) G T (-x), 




B : 


A(x) G T {x) G T {-x) G T {-x) , G T {x) G T {x) 


A(-x) G T (-x 


C : 


A 2 (x) G T (-x) G T (-x) , G T (x) G T (x) A 2 (- 


-x), 


D : 


A(x) G T (x) A(-x) G T (-x) , 




E : 


A 2 (x) A(-x) G T (-x) , A(x) G T (x) A 2 (-x] 




F : 


A 2 (x) A 2 (-x). 





(4.2) 

As the product A(x)A(—x) of zero-temperature propagators involves the combination 
Q(xi)Q(— X4), terms of the classes D,E and F vanish identically. The maximum 
number of B-functions a given term can contain - in order not to be zero - is two. 
Also, the arguments of the two G-functions have to coincide as it is the case with the 
terms of class C. We conclude that we only have to consider the cases A, B and C. 

Contributions from the classes A and B are related to integrals over the torus, 



J d ds+1 x (d r G T d r G T d s G T d s G T + Ad r A d r G T d s G T d s G 7 
which are not singular at d s = 2. 



(4.3) 



Terms of class C, on the other hand, do lead to an ultraviolet-divergent integral. 
Consider, e.g., the term 



d r A{x) d r A{x) d s G T (-x) d s G T {-x) , 



(4.4) 
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where we now have included the derivatives. In the limit d s — >-2, the zero-temperature 
piece d r A(x) amounts to 



x' 



x 4 2 



X 



47X4- 



d r A(x) oc — - exp — . (4.5) 



The Taylor series of the function d s G T (—x), evaluated at the origin, starts with a 
term linear in x, 

d s G T (-x) = d as G T (-x)\ x=0 x a + 0(x 3 ) . (4.6) 

Inserting this term into Eq. (14.41) . in the limit d s -+2, we end up with the following 
contribution in K, 

K oc Jd 2 xdx^-^j e-* 2/2 ^x 2 oc Jdx 4 ^, (4.7) 

which is logarithmically divergent in the ultraviolet. One also readily checks that this 
is the only term that has to be subtracted: The cubic Taylor term in the expansion 
of d s G T , Eq. (l4.6p . leads to a convergent contribution to the integral K. 

In order to isolate the ultraviolet singularity in Eq. (l4.7p . we first cut out a sphere 
S of radius |«S| < /5/2 around the origin and decompose the integral involving the 
contributions of class C as 

d ds+1 x d r Ad r A d s G T d s G T 

= (d ds+1 xd r Ad r Ad s G T d s G T + / d ds+1 x d r Ad r A d s G T d s G T . (4.8) 
Js Jr\s 

In the limit d s — > 2, the integral over the complement T\S of the sphere is not singular. 
The divergence is contained in the integral over the sphere. Here we subtract the 
singular term ( 14. 7p . which leads us to 
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U ds+1 x d r A{x)d r A{x) d s G T {-x)d s G T (-x) 
fd ds+1 x d r A(x)d r A(x) Q ss (x) 
+jd d ° +1 xd r A(x)d r A(x)d as G T (-x)\ x=0 dp s G T (-x)\ x=0 x a x^ (4.9) 

with Q ss (x) defined as 

Q ss (x) = d s G T (-x)d s G T (-x) - d as G T (-x)\ x=0 d(j s G T (-x)\ x=0 x a x? . (4.10) 

The first integral on the right hand side in Eq.f l4.9p is now convergent. The second 
integral, however, is divergent in the ultraviolet. The last step in the isolation of this 
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divergence consists in decomposing the singular integral as follows: 

Jd d ° +1 xd r A(x)d r A{x)d^G T (-x)\ x=0 dp s G T (-x)\ x=0 x a xP 

= [ d d * +l x d r A(x) d r A(x) d as G T (-x)\ x=0 dp s G T {-x)\ x=0 x a x? 

-[ d d ° +1 xd r A(x)d r A(x)d as G T (-x)\ x=0 dp s G T (-x)\ x=0 x a xP . (4.11) 
Jn\s 

Now the ultraviolet singularity is contained in the integral over all Euclidean space, 
which can be written as 

/ d d ° +1 x d r A(x) d r A(x) d as G T (-x) \ x=0 d Ps G T (-x)\ x=Q x a x p 

J K 

= W±j) T d 3 +2 { R) 

nQj it; 3d s ids+i Vr" > 

In the limit d s 2, the above regularized expression is divergent because the T- 
function develops a pole. Finally, we subtract this divergent expression from the 
integral K and define the renormalized integral K as 

K = K-2 fd i xd r A{x)d r A{x)d as G T {-x)\ x=0 d l3s G T {-x)\ x=0 x a x p 

= J d 3 x (d r G T d r G T d s G T d s G T + 4 d r A d r G T d s G T d s G T ^j 

+2 / d 3 x d r A d r A d s G T d s G T + 2 [d 3 x d r A d r A Q ss 
Jt\s Js 

-2 f d 3 xd r Ad r Ad as G T (-x) \ x=0 dp s G T (-x) \ x=0 x a x? . (4.13) 
Jn\s 

While the numerical evaluation of the finite expression K poses no problems and is 
considered in appendix El here we proceed with our discussion regarding the diver- 
gent piece in K. The point is that this divergence originating from the three-loop 
graph 8c can be canceled by an identical divergence in the two-loop graph 8d. More 
precisely, the divergence in the three-loop graph 8c can be absorbed into the next- 
to-leading effective constants l\ and I2 which come along with the two-loop graph 8d. 
This logarithmic renormalization of higher-order effective constants is standard in a 
Lorentz-invariant framework. Indeed, the absorption of ultraviolet divergences occur- 
ring in loop graphs into next-to-leading order effective constants has been discussed a 
long time ago in the context of chiral perturbation theory, i.e. the low-energy effective 
theory of Quantum Chromodynamics [41]. Since the issue is rather technical, we do 
not present it in the main body of the paper - rather we relegate it to appendix [Bl 

Collecting all terms contributing to the free energy density up to order p s , and 
writing the integral K as 

A» = T^, a = ^=^, 7 =£, (4.14) 
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the low-temperature expansion of the free energy density for the two-dimensional ideal 
ferromagnet finally takes the form 



1 _ 9 A e~ an U A e~ an 



- ih T2 Y.— I4 T3 E 

3(4Z| - Cl S 7 ) 



47T7 ^— ' n 2 7rS7 3 ^— ' n 3 

n=l ' n=l 



r 4 E 



7rS 2 7 5 ^— ' n 4 

n=l 

H^Ml^l - - ^0). ,4,5, 



n=l 



The quantities l± and I2 are the renormalized next-to-leading effective constants (see 
appendix [B]) . Note that the next-to-leading effective constant £3, much like the next- 
to-next-to- leading effective constant c\, does not require renormalization - both quan- 
tities are finite as they stand. 



While this series has been investigated before [43M54| . both within modified spin- 
wave theory and Schwinger-Boson mean field theory, we want to emphasize that all 
these references dealt with free magnons. The manifestation of the spin-wave in- 
teraction in the low-temperature expansion for the free energy density of the two- 
dimensional ideal ferromagnet is considered here for the first time to the best of our 
knowledge. Furthermore, our rigorous approach is completely systematic and does 
not resort to any kind of approximations or ad hoc assumptions. 



5 Thermodynamics of Two-Dimensional Ideal Fer- 
romagnets 

Using the representation (I4.15P for the free energy density we now discuss the ther- 
modynamic properties of two-dimensional ideal ferromagnets. Let us first consider 
the low-temperature series for the pressure. Because the system is homogeneous, the 
pressure can be obtained from the temperature-dependent part of the free energy 
density, 

P = z -z. (5.1) 
Accordingly, up to order p 8 , the low-temperature series for the pressure takes the form 

P = 770 T 2 + Vl T 3 + r? 2 T 4 + O(p 10 ) } (5.2) 
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where the coefficients rji are given by 



no 



1 v^e~™ 



E 



4-7T7 ^— ' n 2 

' n=l 

7 OO 



E 



(7(1 



7rE7 3 ' 

' n=l 



3(4/ 2 - C1E7) y> e~ m 

n=l 

8/1 + 6/ 2 + 5/ 



2 



+ 24^4 \E^| +2EV (5 ' 3) 

Note that these coefficients depend on the dimensionless ratio a = fiH/T. Remark- 
ably, the spin-wave interaction starts manifesting itself only at order p 8 oc T 4 through 
the last two terms in the coefficient r]2. While the former contribution, which involves 
the renormalized effective constants l\ and I2, originates from a two-loop graph, the 
latter contribution, proportional to the dimensionless function k(a), comes from a 
three-loop graph. All other contributions in the pressure are related to one-loop 
graphs, i.e. to graphs that describe noninteracting magnons. In the above series for 
the pressure they contribute at order T 2 ,T 3 and T 4 . 

The ratio a = fiH/T in the series for the pressure can take any value, provided 
that the temperature and the magnetic field are both small with respect to the in- 
trinsic scale of the underlying theory. In the present case of the two-dimensional ideal 
ferromagnet, this scale may be identified with the exchange integral J of the Heisen- 
berg model. In the following we will be interested in the limit T ^> [iH which we 
implement by holding T fixed and sending the magnetic field to zero. By keeping the 
fixed temperature small compared to the scale J, we thus never leave the domain of 
validity of the low-temperature effective expansion. 

Formally, the limit a — > poses no problems for the one-loop contributions in the 
pressure. The corresponding coefficients in (I5.3P become temperature independent 
and the sums reduce to Riemann zeta functions. However, one has to be very careful 
by taking this limit in the interaction part contained in the coefficient 7/2. In order to 
address this issue, we first take a closer look at the dimensionless function k(a) which 
is depicted in Fig. [2j 

In the limit a— >-0, this function can be parametrized as 

k{p) = k In a + ki + k 2 a InV + k 3 a In a + k 4 + 0(a 2 InV) , (5.4) 

where the coefficients ki are pure numbers. Here we only need kg and k\ which take 
the values 

k = = 6.818 x KT 4 , ki = 2.18 x 1(T 3 . (5.5) 
4608 
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Figure 2: The function k(a), where a is the dimensionless parameter a = fiH/T. 

The function k(a) thus diverges logarithmically, indicating that the limit a— >0 in the 
above representation for the pressure is rather subtle. In appendix O it is shown that 
the sum of the two- loop and the three- loop contribution in the pressure, 



8li + 6/ 2 + 5/ 3 
24tt 2 S 2 7 4 
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2Y? r 



k{a)T 4 



(5.6) 



remains finite in the limit a—tO and can be written as 



T 4 In 

2S 2 7 3 V T 



24tt 2 E 2 7 4 



C(2)'T 4 . 



(5.7) 



The quantity A p is a scale which depends on the renormalized next-to-leading order 
effective constants l\ and I2, as well as on the coefficients ko and k\. In the absence of a 
magnetic field, the low-temperature expansion for the pressure of the two-dimensional 
ideal ferromagnet hence takes the form 



' C(2)T 2 + 1:1 



47T7 



7rE7 3 



5/, 
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C(3) T 3 + 
k 
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7rS 2 7 5 



C(4) T 4 



T 4 In 
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+ 0(p 



10 \ 



(5. 



Interestingly, we are not dealing with a simple power series. At order p 8 , where the 
spin- wave interaction sets in, a logarithmic contribution T 4 In T emerges. Although 
this term contains a scale related to next-to-leading order effective constants of £ 4 //, 
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the coefficient of this term is fixed by the effective constants S and F of the leading 

'2 

-eff 



order Lagrangian and by the quantity fco which is a pure number. 



In a Lorentz-invariant setting, the occurrence of such logarithmic terms is well- 
known. The low-temperature expansion of the pressure in Quantum Chromodynamics 
involves a term T 8 InT in the chiral limit [§[. In fact, a term T 8 InT also occurs in the 
low-temperature expansion of the pressure of the three-dimensional antiferromagnet 
in the absence of a staggered field Q. Note that antiferromagnetic magnons, much 
like pions in Quantum Chromodynamics, display a linear, i.e. relativistic dispersion 
law. 

While the limit H — > in the interaction contribution turned out to be rather 
subtle, one may have doubts in general about switching off the magnetic field. After 
all, we are dealing with a planar system, where we are facing the Mermin- Wagner 
theorem. However, as shown in Ref. switching off the magnetic field in the low- 
temperature series for the pressure does not really numerically affect the series. We 
will discuss the issue in more detail in the next section where we will consider the 
magnetization. 

Let us now derive the low-temperature series for the energy density u, the entropy 
density s, and the heat capacity cy of the two-dimensional ideal ferro magnet. They 
are readily worked out from the thermodynamic relations 

dP du ds 

s = W , u = Ts-P, <V = - = T-, (5 .9) 



and amount to 
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In the above series, the spin-wave interaction manifests itself in the last two terms 
which originate from two-loop and three-loop graphs. While the former involves the 
renormalized next-to- leading order effective constants l\ and I2, the latter is propor- 
tional to the dimensionless function k(a) and its derivatives. 

If we switch off the magnetic field, these low-temperature series turn into 



5/ 3 
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C(2) 2 T 4 + ^_ t 4 j 3 In (^) - 1 j + 0(p 10 ) , (5.13) 
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(5.14) 



C ^ ^C(2)T + ^C(3)T 2 + 36 ( 4/ 3-f^ C (4)T3 
27T7 7rL7 d 7rL 2 7 tl 

+^lC(2) 2 T' + J-jT' {l21n (£) - 7 } + . (5.15) 

Again, the interaction part does not lead to a simple power of the temperature in the 
limit H — > - rather, it also involves a piece logarithmic in the temperature. 

We want to stress that the structure of the various low-temperature series derived 
in this section is an immediate consequence of the symmetries of the two-dimensional 
ideal ferro magnet. Unlike modified spin- wave theory or Schwinger-Boson mean field 
theory, the effective Lagrangian technique does not resort to any ad hoc assumptions or 
approximations, but is completely systematic. Moreover, to the best of our knowledge, 
the manifestation of the spin-wave interaction in the low-temperature series of the 
partition function of the two-dimensional ideal ferromagnet has never been studied 
before. In particular, to address this problem with spin-wave theory or Schwinger- 
Boson mean field theory up to the order considered in the present analysis, appears 
to be beyond reach. 

At the end of this section, we like to compare the low-temperature expansions for 
the free energy density of the ideal ferromagnet in two and three spatial dimensions, 
respectively. In the absence of an external magnetic field they exhibit the following 
general structure, 

z ds =2 = -^-^-^T^^TMnT + ^p 10 ), 

z ds = 3 = —h Q Ti — h\T^ — h 2 T^ — h 3 T 5 — b^T 11 / 2 + 0(p 12 ) , (5.16) 

where we have highlighted all contributions which are related to the spin-wave inter- 
action. Note that in the case of the two-dimensional ideal ferromagnet, the coefficient 
7)2 has two parts. The former one, r^, contains the free magnon part as well as the in- 
teraction contribution proportional to I3. The latter one, fj^ , is exclusively due to the 
spin-wave interaction which involves the renormalized next-to-leading order effective 
constants l\ and li- 

In three dimensions, the term of order T 5 oc p 10 in the free energy density is the 
famous Dyson interaction term - in the effective framework this is a two-loop effect 
(see Ref. [4|). Corrections to the Dyson term were considered in Ref. [5j, pointing out 
that the leading correction is of order T 11 / 2 oc p u and originates from a three-loop 
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graph. Interestingly, in the case of the three-dimensional ideal ferromagnet, the two- 
loop and the three-loop contributions occur at different orders in the low-temperature 
expansion. As we have outlined in the present article, the same two-loop and three- 
loop graphs, in the case of the two-dimensional ideal ferromagnet, all contribute at 
the same order p 8 and lead to a logarithmic term T 4 In T in the free energy density. 
As argued in Ref. [if, the different organization of Feynman diagrams in three and 
two spatial dimensions is a consequence of the suppression of ferromagnetic loops in 
Feynman graphs: In three spatial dimensions, each loop in a Feynman diagram is 
suppressed by three powers of momentum. In two dimensions, on the other hand, 
ferromagnetic loops are only suppressed by two powers of momentum. 



6 Magnetization and Mermin- Wagner Theorem 

Let us now consider the magnetization. With the expression for the free energy density 
(I4.15p . the low-temperature expansion for the magnetization 

of the two-dimensional ideal ferromagnet takes the form 

^11 = 1 - & T - arT 2 - a 2 T* + <D(p 8 ) . (6.2) 
The coefficients d, depend on the dimensionless ratio a = fiH/T and are given by 



1 ^ e ~ m 



a 



«2 



E 



71=1 

E 



n=l 



7rS 2 7 3 ' n 2 



n=l 

an 



3(4l 2 3 - C1E7) 



7rS 3 7 5 ^— ' n 



n=l 

00 



8Z1 + Qh + 5/ 3 e~™ ^ e-™ 1 dk{a) 

12tt 2 S 3 7 4 ^ m ^ n 2 2S 3 7 3 da ' 1 ' ' 

m=l n=l 

Obviously this series has a problem if we want to switch off the magnetic field, as the 
leading coefficient do then diverges. The point is that it is completely inconsistent to 
take the limit a — > in the low-temperature expansion of the magnetization, as we 
will discuss now. 

We should keep in mind that the present effective calculation is based on the 
assumption that the internal spin symmetry 0(3) is spontaneously broken. While 
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this assumption is fulfilled at zero temperature, at finite temperature, however, spon- 
taneous symmetry breaking in the two-dimensional Heisenberg ferromagnet cannot 



occur according to the Mermin- Wagner theorem |55j. Rather, at finite temperature, 
the low-energy spectrum of the two-dimensional ideal ferromagnet exhibits a nonper- 
turbatively generated energy gap and the correlation length of the magnons, although 



still exponentially large [56 



i 



- C,aS 2 \l — exp 



(6.4) 



T 



no longer is infinite. 



We have mentioned before that these nonperturbative effects are so tiny that 
they cannot manifest themselves in the low-temperature series for the free energy 
density. Likewise, the series for the pressure, the energy density, the entropy density, 
and the heat capacity derived in the previous section are also valid as they stand 
- the subtleties raised by the Mermin- Wagner theorem are not relevant for these 
thermodynamic quantities. 

However, for the magnetization matters are quite different. There, the nonpertur- 
batively generated energy gap does not lead to tiny corrections in the low-temperature 
expansion - rather these effects become the dominant ones. We can easily see this by 
expanding the series in the small parameter a = fiH/T. While the free energy density 
amounts to 

z = -EfxH - ^-T 2 |c(2) + al n( r-a-^ + 0(a 3 ) } + 0(T 3 ) , (6.5) 

the magnetization takes the form 

£(T,tf) = £ + i i-T{lna-! + ^ + 0(a 4 ) } + 0(T 2 ) . (6.6) 

In the effective theory perspective, the quantity AEh, 

AE h = liH, (6.7) 

is the energy gap induced by the weak external magnetic field - this is the quantity 
that appears in the parameter a. However, there is a further mechanism contributing 
to the energy gap, the one that generates a gap nonperturbatively. In analogy to the 
definition of the correlation length, related to the energy gap AEh [H , 



the nonperturbatively generated energy gap AE np is connected to the correlation 
length Cn P by 

_7 

Sn 



AE np = . (6.9) 
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In our low-temperature series, so far only the energy gap AE H is accounted for in the 
parameter a. 



We have to remember that we implement the limit T 3> /iif by holding T fixed 
and sending the magnetic field to zero. Keeping the fixed temperature small compared 
to the underlying scale given by the exchange integral J, we never leave the domain 
of validity of the effective low-temperature expansion. Now for a fixed value of the 
temperature, the nonperturbatively generated energy gap AE np is just a constant. 
However, unlike the energy gap induced by the external magnetic field, the generation 
of AE np is a dynamical effect we cannot control. In particular, while we can switch off 
the magnetic field and thus AEh, we cannot switch off AE np - it would be inconsistent 
to consider values for a smaller than a np = AE np /T in the effective low-temperature 
series. 

At very small temperatures - let us say T/J = 1/100 - the relations (16.41) and 
(16.81) imply that the condition 

AE np = AE H (6.10) 
is satisfied if the ratio fiH/T takes the value 

^ = ^1 e" 100 ^ 10- 131 (£ = §). (6.11) 

This tiny ratio thus leads to a negligible contribution to the free energy density ( 16. 5 p - 
here, although conceptually inconsistent, taking the limit a — > does not numerically 
affect the series. However, in the magnetization (16.61) the effect of the nonperturba- 
tively generated energy gap is rather large because of the logarithm In a - taking the 
limit a — > in the magnetization would be completely inconsistent. 

In conclusion, at nonzero temperature the nonperturbatively generated correla- 
tion length is finite and there is always an energy gap AE np in the spectrum of the 
two-dimensional ideal ferromagnet. The question is whether or not this gap is numer- 
ically relevant in the low-temperature series - for the thermodynamic quantities it is 
irrelevant, for the magnetization, however, it is indeed relevant. 

The low-temperature expansion of the magnetization of two-dimensional ideal fer- 
romagnets has been considered before within the formalism of double-time-tempera 



ture Green functions and spin-wave theory. The explicit expressions given in Ref. [53 



and Ref. [54|, respectively, are consistent with our effective analysis provided that we 



express the effective constants 7 and I3 in terms of microscopic constants as |1| 

l = JSa\ h = ^f- (6.12) 

Here J is the exchange integral of the Heisenberg model and a is the distance between 
the sites on the square lattice. 
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We have to emphasize, however, that the references 53|, |54| were restricted to free 
magnons. The effect of the spin-wave interaction on the low-temperature expansion 
of the magnetization of two-dimensional ideal ferromagnets has been considered for 
the first time in the present work. 



7 Conclusions 



Using the method of effective Lagrangians, we have evaluated the partition func- 
tion of the two-dimensional ideal ferromagnet up to three loops and have derived the 
low-temperature series for various thermodynamic quantities including the magneti- 
zation. In particular, we have shown that in the absence of an external magnetic 
field, the spin-wave interaction starts manifesting itself in the form of a logarithmic 
term T 4 InT in the free energy density. To obtain this result, we had to renormal- 
ize and numerically evaluate a specific three-loop graph which turned out to be the 
piece of resistance. Much like in the case of three-dimensional ideal ferromagnets, the 
spin-wave interaction is also very weak in two-dimensional ideal ferromagnets. 

We have discussed in detail the implications of the Mermin- Wagner theorem for 
the low-temperature series derived in the present work. While the series for the free 
energy density, pressure, internal energy density, entropy density and specific heat 
capacity are also valid if the magnetic field is switched off, one has to be careful in the 
case of the magnetization where the effect of the nonperturbatively generated energy 
gap cannot be neglected. 

Although various authors have considered the low-temperature properties of two- 
dimensional ideal ferromagnets before within spin- wave theory, Schwinger-Boson mean 
field theory and double-time-temperature Green functions, to the best of our knowl- 
edge, they all restricted themselves to free magnons. In particular, none of these 
authors identified logarithmic terms in the low-temperature expansion of the parti- 
tion function. The systematic effective field theory method thus proves to be more 
powerful than conventional condensed matter methods where a three-loop analysis 
appears to be beyond reach. 

Such logarithmic terms are well-known in particle physics, i.e., in chiral perturba- 
tion theory which is the effective field theory of the strong interaction described by 
Quantum Chromodynamics. They also occur in the context of the three-dimensional 
antiferromagnet, whose spin waves obey a linear (relativistic) dispersion law. These 
logarithmic terms are a consequence of the structure of the ultraviolet divergences 
in the Goldstone boson propagator. In Quantum Chromodynamics, in the case of 
three-dimensional antiferromagnets as well as two-dimensional ideal ferromagnets, 
these divergences can be absorbed into next-to-leading order effective constants by 
logarithmic renormalization. Interestingly, in the case of the three-dimensional ideal 
ferromagnet, the ultraviolet divergences of the propagator do not require logarithmic 
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renormalization [5] : The low-temperature series for the thermodynamic quantities and 
the magnetization simply involve integer and half-integer powers of the temperature. 



The present study regarding the low-temperature properties of two-dimensional 
ideal ferromagnets is on the same footing as the analysis performed in Ref . jij , which 
dealt with three-dimensional ideal ferromagnets. In both complete and sys- 

tematic analysis of the partition functions was given up to three-loop order. In the 
three-dimensional case, the effect of the spin-wave interaction on the partition func- 
tion has been considered by numerous authors before, in particular by Dyson in his 
monumental work jij . The effective three- loop analysis presented in Ref. [f| went one 
step beyond Dyson. In the case of the two-dimensional ideal ferromagnet, however, 
to the best of our knowledge, the effect of the spin-wave interaction on the partition 
function has never been discussed so far. The present study attempts to fill this gap 
in the literature. 
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A Numerical Evaluation of the Cateye Graph 

In this appendix, we consider the numerical evaluation of the three-loop graph 8c. 
The various terms in the relevant expression K, 

K = Jd 3 x(d r G T d r G T d s G T d s G T + 4:d r Ad r G T d s G T d s G T ^j 

j3„ a a A 3 fiT q nT , r> / j3. 



+2 / d 6 x d r A d r A dsG 1 dsG 1 + 2 / d 6 x d r A d r A Q ss 
Jt\s Js 

-2 / d 3 xd r Ad r Ad as G T (-x)\ x=0 dp s G T (-x)\ x=0 x a xP, (A.l) 
Jk\s 

only involve the variables r = \x\ and t = x±, such that the integrals become in fact 
two-dimensional, 

d 3 x = 2nrdrdt. (A.2) 

A very welcome consistency check on the numerics is provided by the fact that the 
result has to be independent of the radius of the sphere S. 



It is convenient to introduce the dimensionless integration variables rj and £ 

i It 
2V7 



1 T 

r] = Tx 4 , £=oA/ — Fl- ( A - 3 ) 
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In the integrals over the torus which involve quartic and triple sums - the first two 
terms in Eq. ( lA.ip - we first integrate over all two-dimensional space, ending up with 
one-dimensional integrals in the variable rj. For the quartic sum we obtain 



d 3 x d r G T {x) d T G T {x) d s G T {~x) d s G T {-x) 

+n 2+ n 3+ n 4 ) Q^n^n^) 



32tt 3 7 5 J 1/2 ' ^ 

' ^ V 2 ni... n 4 =l 



-3 



Q(v, ni, n 2 , n 3 , n 4 ) = -7— * -3 , (A.4) 

(77 + ni)(77 + n 2 ){-r) + n 3 )(-r] + n 4 )j 



while for the triple sum we get 



(f^9 r AW <9 r G J (x) dsG 1 (-x) dsG 1 (-x) 

1/2 00 



1 Z" 1 / 2 . . 

—3- T T 4 di? ]T e-^+^)g(77,0,n 2 ,n3,n 4 ) 
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where 



77 ri+n 2 -r?+n 3 -»7+n 4 

(5(77, 0, n 2) n 3 , n 4 ) = y> ^—2 , (A.5) 

^7(77 + n 2 ){-T) + n 3 )(-77 + n 4 )J 

In the case of the triple sums the integration over 77 only extends over the interval 
[0, |], due to the O-function contained in the zero-temperature propagator A(x). 

Note that the quantities Q(r], m, n 2 , n 3 , n 4 ) and 0, n 2 , n.3, n 4 ) depend in a 
rather nontrivial manner on the summation variables. The slowest convergence for 
the expressions Eq. (IA.4[) and Eq. ( 1A.5j) is observed for a = 0, because no exponential 
damping occurs. The numerical summation has been performed in a " Cartesian" way 
as follows. We first define the vector Ni = (ni,n 2 ,n 2 ,ri4). The first partial sum Si in 
the quartic series simply corresponds to the combination Ni = (1, 1, 1, 1) of indices. 
The second partial sum 5*2 then contains all combinations of indices in the vector N 2 
with at least one index equal to two: (2, 1, 1, 1), ... , (2, 2, 2, 2), etc. For large values 
of i and for a — 0, the partial sums Si converge according to 1/Sj 2 . Proceeding in an 
analogous manner for the triple sums, the asymptotic behavior turns out to be the 
same. 
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Expressions suitable for the numerical evaluation of the remaining three integrals 
of Eq. flA.ip involving double sums are 

d 3 xd r A(x) d r A{x) d s G T (-x) d s G T (-x) 

T\S 

1 pS poo 00 

_£2( 2 + 1 , 1 ) 

P(£,r},n 1 ,n 2 ) = ^ , (A. 7) 
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d 3 xd r A(x)d r A(x)Q ss {x) (A.8) 
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(A.9) 
and finally, 
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|77 2 nin 2 | 



In the above integrals the radius of the sphere has been chosen as S = |. For large 
values of i and for a = 0, the partial sums Si related to the above three expressions 
involving double sums converge according to 1 /Si 2 , i.e. the asymptotic behavior is 
the same as for the triple and quartic sums. 
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B Logarithmic Renormalization of the Effective Con- 
stants l\ and I2 



In this appendix we discuss how the ultraviolet divergence in 

/ d d ° +1 x d r A(x) d r A(x) d as G T (-x) \ x=0 dp s G T (-x)\ x=0 x a x? 
Jn 



-,j , j- Ms. Ms+I \l I I / j d s +2 I \ 9 

2 Ms+ °7T 2 >y 2 y n=l n 2 / A 



originating from the three-loop graph 8c, can be absorbed into the next-to-leading 
effective constants l\ and 1% which show up in the two-loop contribution z$r de i, 



T 



24vr 2 S 2 7 4 ' \ ^ " 

' \ n=l 

The singularity in (IB . 1|) is due to the T-function which contains a pole at d s = 2, 

T(l - y ) = -1E + 0(d s - 2) , (B.3) 

where je is Euler's constant. 

Now the two expressions (IB. ip and (1B.2p have the same structure. They both 
involve the same infinite series and are both proportional to four powers of the tem- 
perature. Taking into account the prefactor of the integral K from Eq. (I3.23p . they 
add up to 

The pole in the T-function can thus be absorbed into the combination 8/1 + 6/2 + 5/3 
of next-to-leading order effective constants. Note however, that Z3 does not require 
renormalization. This effective constant already occurred in the one-loop result for 
the free energy density (13.1 j) and is perfectly finite. The combination 8/1 + 6/2, on 
the other hand, is divergent and absorbs the pole in the T-function as we now show 
explicitly. 

Introducing an arbitrary renormalization scale pL, which should not be confused 
with the symbol // denoting the magnetic moment, the divergent quantity A, 

GL rfs-2 

A = T(l - |) (pH)^ , (B.5) 



d 3 -2 



can be decomposed into 

l] .: tl 

(//#) 2 -p,— . (B.6) 
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A = r(i-|)pV + r(i-| 



The first term contains the singularity. Note that it does not involve the magnetic 
field, but depends on the arbitrary renormalization scale ft. The second term is not 
singular in the limit d s — > 2, but approaches the finite value 

-inf^), (B.7) 



such that the sum (IB.4[) takes the form 



87!Gi) + 67 2 (£) + 5ls 3 ^ ,fiH\ j y4 e ~" HnP Y (B 8) 
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The quantities and I 2 (/i) are the renormalized, i.e., finite effective constants 

which depend on the renormalization scale ft. It is important to note, however, that 
the curly bracket is independent of this scale. Moreover, if we chose the renormaliza- 
tion scale as 

fi = fiH=^, (B.9) 

then the logarithm drops out and it is understood that l\(p) and \ {p) are the renor- 
malized effective constants evaluated at this specific choice of the scale. The choice 
that jx be one tenth of the underlying scale given by the exchange integral J, is moti- 
vated by comparison with chiral perturbation theory. There, in order for the logarithm 
to drop out, the scale is to be identified with the pion mass M, which is about one 
order of magnitude smaller than the underlying scale Aqcd- 

From a more practical point of view we can say that the net effect of this whole 
renormalization procedure is simply this: The divergence in (IB. 1[) can be taken care of 
by rewriting the expression z^j de ] in terms of the two renormalized effective constants 

li and I2 evaluated at the renormalization scale /2 = fiH as 



T 

8[<fe] 



8Z1 + 6Z 2 + 5/3 ^ 4 I 



x n=l 



C Pressure at Zero Magnetic Field 

In this final appendix, we discuss the low-temperature expansion for the pressure at 
zero magnetic field. While this limit in the one-loop contributions formally poses no 
problems, switching off the magnetic field in the interaction contribution of order p 8 , 

2 



8h + 6Z 2 + 5Z 3 e~™ ^ 1 ^ LiH 
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is rather subtle. Using the representation (15. 4p for the dimensionless function k(a), 
the interaction contribution - in the limit H — y - amounts to 



k J 7T 2 (8/1 + 6/2 + 5/3) +ln fv£\ + h 1 T 4 / C2 x 



2E 2 7 3 1 432 7 fc IT/ fc 

The effective constants /1, 12 and Z3 have the same dimension as 7 = F 2 /£. For reasons 
which become obvious in a moment, we write the renormalized effective constants as 



7 , / Ai \ is - ( A 2 \ 12 432fe 

/1 = 7 c In — , l 2 = 7 c In — , c = — — 



7T 2 



(C.3) 



where we have introduced the two scales Ai and A2 which are independent of the 
magnetic field. Inserting this representation for l\ and I2 into the first term in the 
curly bracket of the expression (lC.2p . one obtains 

M ia ^ +in ^ + ^} T '' {ca) 

where the scale A p is given by 

A p = v^As" e fel / fc0 w 24.5 y/KJT 2 . (C.5) 

The essential point is that in the sum (IC4[) the magnetic field drops out. We conclude 
that, if one switches off the magnetic field in the low-temperature expansion for the 
pressure, the interaction contribution does not just manifest itself in the form of a 
simple power of the temperature - rather, it also involves a logarithmic contribution 
which reads 

k ° TUn(^). (C.6) 



2S 2 7 3 V T 
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